Interference Alignment for Partially 
Connected MIMO Cellular Networks 



Liangzhong Ruan, Student Member, IEEE, 
Vincent K.N. Lau, Fellow, IEEE, Xiongbin Rao 



Abstract 

In this paper, we propose an iterative interference alignment (lA) algorithm for MIMO cellular 
networks with partial connectivity, which is induced by heterogeneous path losses and spatial 
correlation. Such systems impose several key technical challenges in the lA algorithm design, 
namely the overlapping between the direct and interfering links due to the MIMO cellular topology 
as well as how to exploit the partial connectivity. We shall address these challenges and propose 
a three stage lA algorithm. As illustration, we analyze the achievable degree of freedom (DoF) of 
the proposed algorithm for a symmetric partially connected MIMO cellular network. We show that 
there is significant DoF gain compared with conventional lA algorithms due to partial connectivity. 
The derived DoF bound is also backward compatible with that achieved on fully connected K-pair 
MIMO interference channels. 

I. Introduction 

Recently, there are intense research interests in the area of interference channels and 
the associated interference mitigation techniques. In particular, lA approach can achieve 
the optimal degree of freedom (DoF) in K-pair interference channels [IJ as well as 2-pair 
MIMO-X channels [2]. In |l3l, the lA approach is extended to cellular OFDMA systems 
by exploiting some problem-specific structure such as the channel states being full-rank 
diagonal matrices. In lIH, [|5]|, lIH, [|7]|, the authors extend the lA approach to MIMO cellular 
networks. However, these works have focused on two-cell configuration with one data stream 
for each mobile (MS) JH, (Si or with no more than two MSs in each cell Q, |l7l, and their 
extension to general MIMO cellular networks (with arbitrarily number of cells, MSs and data 
streams) is highly non-trivial. Furthermore, in all these works, a fully connected interference 
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topology is assumed. In practice, we might have heterogeneous path losses between base 
stations (BSs) and MSs as well as spatial correlation in the MIMO channels. These physical 
effects induce a partially connected interference topology. Intuitively, partial connectivity 
in interference topology may contribute to limiting the aggregate interference and this may 
translate into throughput gains in interference-limited systems. Yet, in order to exploit this 
potential advantage, it is very important to incorporate the partial connectivity topology in the 
lA algorithm design. In this paper, we are interested to study the potential benefit of partially 
connectivity in MIMO cellular networks with general configurations and quasi-static fading. 
There are some key technical challenges that have to be addressed. 

• Challenges inherent to MIMO cellular networks: The existing iterative lA algorithm 
designed for interference channels fSl exploit the statistical independency of the direct 
links and the cross links. However, for MIMO cellular networks, there is overlapping 
between the direct links and the cross links as illustrated in Fig. [T] As a result, brute 
force application of the conventional lA schemes in MIMO cellular systems may not 
have desirable performance. 

• Challenges to exploit Partial Connectivity: In practice, MIMO cellular systems are 
usually partially connected due to path losses and spatial correlation, as illustrated in 
Figure [IB. Designing an lA algorithm which can exploit the benefit of partial connectiv- 
ity in the general case is highly non-trivial. While part of this issue has been addressed 
in our prior work BU, the algorithm proposed in |[9l cannot be directly extended to the 
cellular case due to the specific challenges induced by the cellular typology. 

• Challenges due to Quasi-Static Fading: For quasi-static interference networks, the lA 
design may be infeasible [[TOll . However, the lA feasibility checking algorithm proposed 
in [fTOll involves huge complexity of 0{2^^), where is the total number of nodes 
in the network. Such a complexity is intolerable in practice. Hence, a low complexity 
algorithm for checking the lA feasibility conditions on a real-time basis is needed. 

In this paper, we will tackle the above challenges by proposing a novel lA algorithm 
that exploits the partial connectivity topology in MIMO cellular networks. We adopt an 
optimization-based approach and decompose the problem into three sub-problems which 
allows us to tackle the challenges due to MIMO cellular topology and the partial connectivity 
separately. Moreover, we propose a low complexity lA feasibility checking algorithm that 
has worst case complexity of 0{N'^) only. Based on the proposed scheme, we derive an 
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achievable bound of the DoF in a symmetric partially connected MEMO cellular network. We 
show that using the proposed algorithm, the partial connectivity can be exploited to increase 
the total DoF in the MIMO cellular networks. Finally, the proposed scheme is compared with 
various conventional baseline algorithms via simulations and is shown to achieve significant 
throughput gain. 

The following notations are used in this paper: a, a. A, and A represent scaler, vector, 
matrix, set/space, respectively, in particular, R, C represent the set of real number and 
complex number, respectively. The operators (■)^, (■)^, rank(-), trace(-), | ■ |, and dim(-) 
denote transpose, hermitian, rank, trace, cardinality (of a set) and dimension (of a space), 
respectively, span ({a}) denotes the linear space spanned by the vectors in {a}. span({A}) 
represents the space spanned by the column vectors of A. Q{S, N) denotes the Grassmannian 
IfTTTl . which represents the space of all the S dimensional subspaces of the dimensional 
space. 

IL System Model and Problem Formulation 

A. MIMO Cellular Networks 

We consider a MIMO cellular system with G BSs, each of which serves K MSs, as 
illustrated in Fig. [1^. Denote A'^, N^j^ as the number of antennas at BS-(^ and the A;-th MS 
of BS-(7, respectively. Denote dg^ as the number of data streams transmitted to the k-\h MS 
from BS-(/. The received signal at the A;-th MS of BS-f? is given by: 

ygk = \J%\y2^gk,nY,^ni^ni + Z J , Vfc G { 1 , . . . , (1) 

\n=l i=l / 

where Hgfc,„ G C^s*^^^" is the channel state information (CSI) from BS-n to the k-ih. MS 
of BS-^f, x„j e C"''"^\ V„i e C^''^"'"' and G £N^^^d„i information symbols, 

the precoding matrix and the decorrelator matrix, respectively, for the z-th MS of the BS-(7. 
z G <C^ak^^ is the white Gaussian noise with unity variance. The CSI matrices {Ylgk,n\ 
are assume to be quasi-static and mutually independent random matrices. Furthermore, we 
normalize the precoding matrix and the transmit symbols as trace (V^V„j) = dni and 



E 



^f^i trace(x^jX„j) = P„ so that the total transmit power from BS-n is P„. 



B. Partial Connectivity in MIMO Cellular Networks 

We first describe the statistical model of the CSI matrices {Hg^^}. This model was first 
proposed in [121 and widely adopted in literature. 
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Assumption 2.1 (Channel Fading with Dual-sides Correlation): We consider a channel model 
that incorporates both transmit and receive spatial correlation and channel gain, thus: 

Hgfc n = G'gfc.riA^ Hg^ B„fc „ (2) 



'-gk,n — '^gk,n-^gk,n^gk,n^gk,n 

gk,n 



where H^^ „ G C^^'^''^" contains i.i.d. CA/'(0, 1) entries, Gg^^n G M+U{0} is the square root of 
channel gain, Agfc „ G C^gk^'^gk^ ^gk,n ^ C^"^^" represent the receive and transmit spatial 
correlation, respectively. Here A^^Agfc„, B^^B^^^^ are positive semi-definite matrices, 

||AgA:,n||F = ||Bgfc^„||i;' = 1. ■ 

Based on the statistical model of the CSI matrices, we formally define the notion of partial 
connectivity below. 

Definition 2.1 (Partial Connectivity): We define the partial connectivity between BS-n and 
the A;-th MS of BS-g to be the null space of the spatial correlation matrices Ag^ „, B^fc „ 
time channel gain factor Ggk^n- 

• Transmit partial connectivity: A/'^;. „ = ^{Ggk^n^gk,^)- 

• Receive partial connectivity: A/'J^, „ = M{Ggk,n-^gk,n)- ■ 
Remark 2.1 (Physical Meaning of Partial Connectivity): The partial connectivity actually 

describes the effective subspaces of the channel matrices between BSs and MSs in the 
network. For instance, {J^gk -^gk n} represent the subspaces that cannot be perceived by 
the BSs and the MSs, respectively. Hence, the partial connectivity topology of the MIMO 
cellular network is parameterized by the null spaces {J^gk,ni-^gk,n}- Also note that both the 
inter-cell links (i.e. g n) and the intra-cell links (i.e. g = n) may be partially connected. ■ 
We consider a few examples below (as shown in Fig. [l^) to illustrate how the partial 
connectivity model in Definition 12.11 corresponds to various physical situations. Note that 
CSI matrices ligk,n G C^^^ are modeled by ©. 

• Fully connected MIMO cellular network: If Ggk^n ^ and Agk^n, ^gk,n are full rank, 
we have J^g^^n = ■^gk,n = {0}. V(7, k and this corresponds to a fully connected network. 

• MIMO cellular network with spatial correlation: As an illustration, H2i,i has spa- 
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tial correlation such that A2i,i = 

span([0,in,Ar2Vi = span([l,on. 
MIMO cellular network with heterogeneous path losses: Suppose the path loss from 
BS-1 to the second MS of BS-2 is 60 dB and the transmit SNR is 40 dB. Since the 
interference power from BS-1 is negligible compared with the gaussian noise, we can 
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effectively assume 6*22,1 = 0, which gives A/'22,i = J^22,i = illustrated in Fig. [T^. 

C. Stream Assignment and Transceiver Design under Interference Alignment Constraints 

We assume all the BSs in the MIMO cellular network share global CSI knowledge!^ 
{Hgfc „}. We adopt the lA approach to maximize the network total DoF, which is defined 
by D = limsNR-5>oo log(SNR) ' where C is the network sum throughput and SNR is the 
signal to noise ratio. Note that C = D log(SNR)+c'(log(SNR)), DoF gives a first order 
estimation on network throughput. Moreover, it offers some first order simplification to the 
complex throughput optimization on MIMO interference network. Specifically, we would 
like to jointly optimize the data stream assignment {dnj}, precoders {V„j} and decorrelators 
{Unj}, n E {1, G},j E {1, K} policies to maximize the total number of data streams 
Yln=i ^k=i ^nj under the lA constraints^, i.e.: 

Problem 2.1 (lA for MIMO Cellular Networks): 

G K 

max > > dnj (3) 

{d„,},{V„,},{U,fe}^^ 
n — 1 J — i 

s.t.: rank(Ufi,Hgk,gVgk) = dgk, (4) 
U^Hgfc_„V„j = 0, (5) 
trace (V^^.V„j) = dnj, (6) 
E {0,1,...,^^^^^}, Wg,nE{l,...,G}, k,j E {l,...,K}, {n,j) ^ {g,k) 

where d^^^ is the maximum number of data streams for the concerned MS. Constraint ([H) 
ensures that all the direct links have sufficient rank to receive the desired signals while 
constraint ([5]) ensures that all the undesired signals are aligned. 

III. lA FOR Fully Connected MIMO Cellular Networks 

In this section, we shall first solve Problem 12.11 for fully connected MIMO cellular net- 
works, i.e. dim(A/];-, „) = dim(Ar;fc „) = 0. 

'Global CSI is easy to obtain when the network size is small. When the networks size is large, the paitial connectivity 
can be exploited to achieve scalable CSI feedback schemes. For instance, by utilizing heterogeneous path loss, in 1131 , the 
authors propose a scalable CSI feedback scheme for MIMO cellular networks. 

^Under the lA constraints (|5]l, the number of data streams equals to the DoF of the network. 
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A. The Unique Challenge for MIMO Cellular Networks 

In the literature, a common approach towards lA for interference channel is based on the 
interference leakage minimization iteration |[8l. While this approach is designed for standard 



lular network as below: 
Alternatively update pre- 



I 



interference channels, one can extend the framework to MIMO ce 
Algorithm 1 (Extension of Existing Iterative lA Algorithm /lH/j; 
coders V„j and decorrelators IJgk by minimizing the total interference leakage expressions 
in d?]) and ([8]) until the algorithm converges. 

G K 

mill Yl E trace ((Uf,H,,,„V„,)^(Uf,H,,,„V„,)) (7) 

v„^-ec^" g=l k=l 

V-ff.V„.=I (9,fc)#(n,i) 

G K 

^j.^.. E E trace ((Uf,H,,,„V„,)^(Uf,H,,,„V„,)) (8) 

Vn,^7e{l,...,G},j,A;G{l,...,iC}. 

Fig. [5] illustrates the performance of the naive algorithm in a 3-BS fully connected MIMO 
cellular network with K = 2, Nj^ = 5, N^^ = 2, = 1, G {1, 3}, A; G {1, 2}. It is shown 
that the naive algorithm could achieve a total DoF of 3, which is only half of the achievable 
DoF lower bound given in ifTOl . which demonstrates that naive extension of standard iterative 
lA algorithm can perform poorly in MIMO cellular networks. This problem is due to the 
direct link - cross link overlapping issue defined below: 

Definition 3.1 (Direct Link - Cross Link Overlapping): In an interference network, denote 
the set of the channels that carry the desired signals and undesired signals as and H*^, 
respectively. If n ^ 0, then the network has direct link - cross link overlapping. ■ 

As illustrated in Fig. [iK, in conventional MIMO interference network, = {Hmm} 
and H'" = {Hmn '■ 7^ n}, where m, n are the indexes for transmitters and receivers, 
respectively. Obviously, in this case, H^flEI'" = and there is no overlapping issue. However, 
as illustrated in Fig. [l^, in MIMO cellular network, when the number of MSs per cell 
K > \, the intra-cell links also carries over undesired signals and hence = {H^^ g} and 
= {Hgfc,n}- In this scenario, we have that nif'' = 0. As the channel states in 

'The algorithm proposed in (8) is important as it offers a systematic way to obtain lA transceiver design for MIMO 
interference networks with general configuration. Most other existing lA algorithms applies to MIMO interference networks 
with simple specific configuration only. Please refer to (9l Sec. I] for details. 
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appear in ^ and ([8]), when we update the precoders and decorrelators via (|7]) and (O, 
we may also reduce the dimension of the signal space for the desired signals as well. 

B. Problem Decomposition 

We decompose the original problem, i.e. Problem [ZTI into the following three subproblems: 

Problem 3.1 (Stream Assignment): 

G K 

max ^ ^ dnj (9) 



n=l j=l 



K 



S.t. dgkdn,< J2 dn,{K-J2dnk)+ dgk{N^g, - dgk) (10) 

(9,fc)gSjj,(n,j)eS^ (nj)gg^ k=l {g,k)&u 

VSy, §t7 C §, where S = {{g,k) : g e {1, ...,G},k E {1, K}} 
Problem 3.2 (Inter-cell Interference Suppression): 

G K 

min ^ t^^^^ ((Uf,H,,,„V^^.)''(Uf,H,,,„V^^.)) d D 

S.t. V^^. = v:^^. + S„V„^^, V^^. 6 C(<-S-^'^-)x<. (12) 
Uf,U,fe = I,U,feGC^«^^x"'= (13) 

V(/ G {1, G}, A; G where: {d'^j} are the solutions of Problem 13. 1[ matrices 

G C^"X'^"J, Sn G C^''x^^''^^*^=i are isometry matrices whose row vectors combined 
together form a basis for C^"^^, i.e. 

[V^i,V^2,...V^^,S„]''[V^i,V^2,.-.V^^,S„] =1. (14) 

Problem 3.3 (Intra-cell Interference Suppression): 

K 

min Y trace (((U:,)^H„,,„V„,)^((U:,)''H„,,„V„,)) (15) 

S.t. (V„,)^V„, = 1 (16) 
rank ((U^O'^H^.^.V^^.) = rf;^. (17) 

span([V„i,V„2,...,V„;,]) C span ( [Vf,, V^^, V^*^] ) (18) 

where {V^*} and {U*j}, n G {l,...,^},^ G {1,...,/^}, are the solutions of Problem [XIl 
span(X) denotes the linear space spanned by the row vectors of X. ■ 
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Theorem 3.1 (Equivalence between the Original Problem and the Subproblems): For fully 
connected MIMO cellular networks with i.i.d. channel matrices {H^^n}, the optimizing 
variables of Problem [IH is given by {{d*^}, {U*;,}, {V*^}) with probability 1, where {d*^}, 
{U*^.}, {V*^.} are the solutions of Subproblem I3.1[ \3?2\ 13.31 respectively. Furthermore, the 
optimal value of Problem 2.1 is D* = Ylf=i ^nj- 

Proof: Please refer to AppendixEA] for the proof. ■ 

Remark 3.1 (Roles of the Three Subproblems): 

• Problem l3.il determines the stream assignment {dnj} to maximize the sum of the data 
stream numbers (i.e. DoF) of the network, conditioned on the network being lA feasible. 

• Problem 13.21 updates the intermediate precoders {V^^} and decorrelators {Ugk} to 
suppress the inter-cell interferences. 

• Problem \3.3\ further adjusts the precoders {V„j} to suppress the intra-cell interferences. 
Note that after separating the process of inter-cell and intra-cell interference mitigation in 

Problem 13.21 and Problem 13. 3[ only inter-cell channel states {Hg^ .„}, g ^ n are involved 
in Problem 13.21 This property is very important to overcome the cross link - direct link 
overlapping issue. ■ 

Remark 3.2 (The Structure of the Intermediate Precoder): Unlike the existing iterative lA 
algorithm, we have introduced an auxiliary variable, namely the intermediate precoder vari- 
ables {Vrij}. From (fT2l) . Yj^j consists of the core space V^^, the free space S„ and the free 
elements V^^ as illustrated in Fig. |2l This precoder structure in the auxiliary variable enables 
us to separate inter-cell and intra-cell interference suppression. ■ 

Remark 3.3 (The Physical Meaning of Equation (fTSi) )." Constraint (fTSi) is introduced to make 
sure that the desirable inter-cell interference alignment property obtained in Problem 13.21 
is still maintained during the precoder updates {V„j} in Problem 13.31 This is because of 
the following. Suppose {U*;,} and {V^*} constitute the solution of Problem 13.21 We have 



(U;,)^H,fc,„ [V^i,V^;,...,V^*^] =0,\/gj^ne {!,..., G}. From m, there must exist a 
matrix R„ G C^" i -^--E" i such that [V„i, V„2, V„^] = [V^*i, V^*, V^*^] R„, 



(U;,)^H,fc,„ [V„i, V„2, ..^Yuk] = iU;kf^gk,n [V^i, V^*2, V^y R„ = ■ R„ = 0.(19) 



Hence, equation (fT9l ) shows that the inter-cell interference alignment property is preserved 



which leads to the following equation 




March 22, 2012 



DRAFT 



8 



C. Solution to Problem \3.1\ ( Stream Assignment Problem) 

Problem 13.11 is a combinatorial problem whose optimal solution {(i*^.} often involves 
exhaustive search with exponential complexity w.r.t. to the total number of MSs GK. For 
low complexity consideration, we propose the following greedy-based solution. 

Algorithm 2 (Greedy Stream Assignment): 

• Step 1 Initialization: Initialize the stream assignment policy to be the number of streams 
requested by each MSs, i.e. dgk = dfi^"", \/g e {1, G}, k e {1, K}. 

• Step 2 Low complexity lA feasibility checking: 

- Denote t'^^j, fg^, n,g G {1, G}, j, A; G {1, K} as the number of the freedoms, 
i.e. free variables in precoder V.^^ and decorrelator Ugk, respectively. Note that the 
number of freedoms in V^^ are given by the number of elements in V^^ and that 
in IJgk are given by the dimension of Grassmannian Q{dgk, N^^), we have 



K 

t 



dnAK - E ^"'^)' ^.''^ = ^ak{Nl, - dgk). (20) 



fc=l 



- Denote Cgk^nj, n,g & {l,...,^}, j G k G as the number of 
constraints required to eliminate the interference from V„j to \Jgk- Set 

Cgk,nj = dnjdgk, \i g^n] Cgk,nj = 0, otherwise. (21) 

- Use the low complexity lA feasibility checking algorithm proposed in AppendixEB 
to check if the system is lA feasible. If the network is not lA feasible, go to Step 
3. Otherwise, let d*^ = dgk, G {1, G}, A; G {1, A'} and exit the algorithm. 

Step 3 : Update dg/^' = dgi^' — 1 and go back to Step 2, where (g' , k') is given by 

/ G K 

{g', k') = arg max ( ^ ^{cgk,nj + Cnj,gk - dgk^nj - c'nj^gk) - {vlk + vlk - v'lk - v'lk) 

n=l j=l 



arg max 



( G K \ 

2E E dn,-{Nl + Nl,-Adgk + 2) 

\ n=l i=l / 

\ (n,j)#(9,fe) / 



(22) 



where {v'^k^v'^gk}^ Wgk,nj^^'nj,gk} dcnotc the number of freedoms and constraints 

given by (|20l) and (1211) . respectively, with d'^f^ = dgk = 1. 
Theorem 3.2 (Property of the Low Complexity lA Feasibility Checking): The lA feasibil- 
ity constraint in (flOl) is satisfied if and only if it can satisfy the low complexity lA feasibility 
checking in Appendi:?iPBl Moreover, the worst case complexity of the proposed checking 
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scheme is 0{G^K^), which is substantially lower compared with the complexity 0{2'^^^^) 
in conventional lA feasibility checking IfTOl . 

Proof: Please refer to lfT4]| for the proof. ■ 

D. Solution to Problem 13.21 (Inter-cell Interference Suppression Problem) 



The following algorithm solves Problem 13.21 by alternatively updating the intermediate 
precoders {V^^} and the decorrelators {U^fc} to minimize the inter-cell interference, i.e.: 

G K 

mm E t^^'^e ((U5H,,,„V^^.)''(Uf,H,,,„V^^.)) , S.t.: equation CH), (23) 



"J 9 = 1 fc=l 

G K 



min E Yl t^^^e ((Uf,H,,,„V^)^(Uf,H,,,„V^)) , S.t.: equation m. (24) 

9-1 k = l 

Algorithm 3 (Alternative Inter-cell Interference Suppression): 

• Step 1 Initialization : Randomly generate VriG {1, ...,G}, j G {1, ...,K}. 

• Step 2 Minimize interference leakage at the receiver side: At the k-th MS of BS-g 

, update Vgk. Ugk{d) = Ud Yln=l,^gYlf=l Pnji'iigk,n'Vij)i'iigk,n'yij)^ , whcrC Ugk{d) 

is the d-th column of IJgk, i^d[A] is the eigenvector corresponding to the ci-th smallest 
eigenvalue of A, d G {1, dgk}. 

• Step 3 Minimize interference leakage at the transmitter side: At BS-n, update 
V^^, J G {1,...,K}: V^^. = -(S^Q„,S„)-iS^Q„,V^^. ,where Q„, = EU,^nEti 

. Repeat Step 2 and 3 until V^^ and U^fc converges. Set V^* = V^^ + S„V^j and 
u;. = U,,. 

Theorem 3.3 (Convergence of Algorithm\3): For fully connected MIMO cellular network 
with i.i.d. channel matrices {Hg^ „,}, Algorithm [3] converges to a local optimal solution of 
Problem 13.21 Note that global optimality is not assured. 

Proof: Please refer to Appendi)<rC] for the proof. ■ 

Theorem 3.4 (Property of {V^*} and {U*^ }).- For fully connected MIMO cellular net- 
work with i.i.d. channel matrices {Hgk^n}, the converged solution of Algorithm [3] {V^*}, 
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U;^., n e {1, G}, J e {1, K}, satisfy 



rank 



(U*i)^H„l,n 



\ L (U*^)^H„A',n 

almost surely. 

Proof: Please refer to AppendixED] for the proof. 



K 



Y,dnj,^ne{l,..,G} (25) 



E. Solution to Problem \3.3\ (Intra-cell Interference Suppression Problem) 



q+l, VV VV g 



We solve Problem 13.31 by the following constructive algorithm. 
Algorithm 4 (Intra-cell Zero-Forcing): Denote = Wf , W^^^, 
where = (U*y)^H„q_„, q G {1,...,A'}. Each BS does the following for every q G 
{1, A'} to calculate the precoders: 

. Step 1: Perform LQ decomposition for [V^* , V^2) ^uk] = L„(g)Qn(?), where 
Qn{q) is an Y^f^^ dnj x Y^f^^ d^j unitary matrix, and L„(q) is a Ylf=i dnj x Xljli dnj 
lower triangular matrix. 
. Step 2: Set V'„^ = [V^* , V^;, V^}^] Q'^{q), where Q'Jg) is a matrix aggregated by 

the last dng columns of Q^{q). 
• Step 3: Perform singular value decomposition for V^^ = A„qS„gB^^, where S„g is a 
Ng X (i„g matrix, A„g and B„q are A'^ x Ng and (i„q x dnq matrices, respectively. Set 



nq 



A' , where A' is a matrix aggregated by the first dnq columns of A 



'-nq- 



Theorem 3.5 (Optimality of {V^j}): For fully connected MIMO cellular network with 
i.i.d. channel matrices {Hgk^n}, the output of Algorithm |4] {V*y}, n G {l,...,^}, j G 



{1, A'}, is the optimal solution for Problem 13.31 almost surely (with optimal value (intra- 
cell interference power) =0). 

Proof: Please refer to Appendi^jPEl for the proof. ■ 

IV. lA FOR MIMO Cellular Networks with Partial Connectivity 

A. Space Restriction on Transceivers 

In the prior work [9], we have shown that in contrast to the conventional cases, partial 
connectivity can be beneficial to system performance in MIMO interference networks as 
it gives us an extra dimension of freedom, namely the interference nulling to eliminate 



H 
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interference!. In particular, we have found by restricting transceivers to lower dimensional 
subspaces in partially connected MIMO interference network, we can eliminate many lA 
constraints at a cost of only a few freedoms in transceiver design and hence extend the 
lA feasibility region. We exploit the idea of subspace constraint to extend the approach in 
Section [ni] to exploit the partial connectivity in MIMO cellular networks. Specifically, we 
impose the following structure on transceivers: 

Definition 4.1 (Transceiver Structure to Exploit Partial Connectivity): 

• Intermediate precoder with dynamic free space: V^^ = V^^ + S^^ V^^, 

• Decorrelator with dynamic linear filter: IJ nj = S^^U^^, 

where G C^"^^"^ V^^. e C^n.>«^-, S^^- G C^n.x('^-+^n,), u^^. ^ C^^-^^"^')^'^-, 
5*^. G {0, 1, ...K - Eti dn,}, G {0, 1, ...iV^ - dn,}. U 

Remark 4.1 (Space Restriction via New Transceiver Structures): Note that span(V;^j) C 
span(V^j) + span(S*y), and span(U„j) C span(S5^j), space restriction is imposed on V^^ 
and U„j by the new transceiver structure. As a special case, when 5"*^- = — J2k=i ^nj, 
S^j = N^j — dnj, the transceiver structure is reduced to that proposed in Section Hill ■ 

B. Problem Decomposition 

Similar to Section [nil the original Problem [2. II is decomposed into three subproblems. The 
data stream assignment subproblem is modified as below. 
Problem 4.1 (Stream Assignment and Subspaces Design): 



G K 



max 



(26) 



S.t. ^ min {dgk, dim ((span(S^fe) n (ATJ^^^)^)) min (dnj, 

dim((span(VS) + span(S^))n(Ar;,J^)) < ^ d^.S^^ + dgkS^fll) 

VSy, §t/ C §, where E = {{g,k) : g e {1, ...,G},k E {1, K}} 

Ki, V^2, V^^, S^] [V^i, V^2, -.vs., S* J = I, V^^. G (AT^. J\ (28) 

(S^^O^'S;,. = I, S:;,. G (AC,-J^, Vn G {1, ...,G},j G {1, ...,K}. (29) 

''please refer to (9] Section III-A] for detailed elaboration on the concept of interference nulling and the difticulties to 
integrate it into lA processing. We shall omit the details here due to page limitation. 
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The second and third subproblems are similar to Problem 13.21 and 13.31 except replacing 
Y^j, S„, and U„j with Vg, S^* , and S;*U^^. respectively in Problem [3j and [331 where 
{V^;}, {S^* }, and {S^*} are the solutions of Problem gU 

Theorem 4.1: (Connetction of the Original Problem and the Subproblems in Partially Con- 
nected Networks) For partially connected MIMO cellular networks, we have with probability 
1, the solutions of Subproblem 14. 1 [ [XIl 13. 3 [ i.e. {rf*^.}, {U*;,}, {V*^}, are also valid solution 
of Problem 12.11 Hence, the performance of the decomposed problems, i.e. Yl^=i XlfcLi '^gfe 
gives a lower bound of that of the original problem. 

Proof: Please refer to Appendi)iFFl ■ 

C. Solution to Problem \4.1\ ( Stream Assignment and Subspaces Design Problem) 

We extend the greedy -based Algorithm [2] to cater for the partial connectivity in Problem |4~T1 
Algorithm 5 ( Greedy-Based Solution for Problem \4.1\i : 

• Step 1 Initialization: Initialize the number of streams as: dnj = min(rank(H„j d^^^^), 
Wne{l,...,G},je{l,...,K}. 

• Step 2 Calculate the common null spaces: At each BS n G {1, G}, calculate the 
intersection of the null spaces of the inter-cell cross links, i.e. A/'„(M) = r](^g^k)eM-^gk 

M C {{g, k): g^nE{l, G}, k E {1, K}}, as follows: 

- Denote M„ = {{g, k) : ^ 0}. Initialize Mn{^) = C<, Mn{{{g, k)}) = A(;,^„, 
and set the cardinality parameter C = 2. 

- For every M C M„ with |M| = C, if all the subsets of M with cardinality (C - 1) 
are not {0}, calculate Ar„(M) = Ar„(M\{(^', A;')}) n A/'({(c/', A;')}), where {g',k') 
is an arbitrary element in Ksub- Update C = G + 1. Repeat this process until 
Af{M) = {0}, VM C M„ with |M| = C or C = |M„|. 

- For every M C M„ with 7V;(M) ^ {0}, set A/;,(M U ({1, K}\M„,)) = A/'„(M). 
At each MS Mgk, calculate A^J^lM') = nneM'^fgk,n^ M C {n : g e {l,...,^}} 
using a similar process. 

. Step 3 Design V,^^- (i.e. span(V^^.)) : At BS n, n G {1, G}, design Y^j, j G {1, ...,K} 
one by one as follows: For the j-th MS of BS-n, 

- Update the number of streams assigned to the j-th MS of BS-n if there is not enough 
signal dimension left, i.e. update dnj = min (^dnj, Ng — dim ((+A:<jV^fc) + -^(Hnj))) ; 

- Design Y^j based on the principles that A) Y^j is orthogonal to the previous 
designed core spaces and is contained by the effective subspace of the direct link. 
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i.e. V^^- C {{+k<j^nk) +-^(Hnj)) ; B) A subspace which belongs to a null space 
J\f{M) with larger "weight" (i.e. Wn{M{M)), defined below) is selected with higher 
priority. 

Wr,iAf{M))= J2 min(rf,fc,rank(H,fc,„)) (30) 

(g,fe)eM 

From the left hand side of (l27l) . this weight is the maximum number of I A constraints 
that one can mitigate by selecting a one dimensional subspace in A/'(M). 
. Step 4 Design §^ and (i.e. span(S^)> span(S^,)): 
At BS n,ne {1, G}, design {Sj,^.}: 

- A. Generate a series of potential §>l^j{d),d E {0, 1, A^* — J2k=i'^rik} with 
dim(§^(d)) = d based on the principles that A) §^ C {{+ke{i,...,K}^nk))^ , B) 
Same as the principle B in Step 3. 

- B. Choose the best §*„j.- Set = Sl,j{d*), where 

G K 

d* = SiTgmax (^dnjd - ^ ^ min (dg^, rank(Hgfe .„)) 

3=1, k=l 

X min |(V^^. + §*„/d)) n iAfi,,y\) ) . (31) 

Similarly, at each MS Mgk, generate SJ^^((i), d G {0,1,..., — dgk} based on principle 
B. Set = El.id*), where 

G K 

d* = argmax (dgkd - ^ ^min {dgk, {Sl^n {J\f^i,^J^\) 

n=l,f^g i=l 

X min |(V^,. + §^(d)) n {X'gk,n)^\) ) . (32) 

• Step 5 lA Feasibility checking: Similar to Step 3 in Algorithm |2l set v^^j = dnjS^^j, 

^gfe = dgkS'g^, where S'*^- and S'^g^ are defined in Definition l4.1[ Set Cgk^nj = min {dgk-, |span(S 
n(A/;V)^l)"^inKi'l(span(V^^0 + span(S^^0)n(A4*,^^^ if g ^ n- Cgk,nj = 0, 
otherwise. Use the low complexity algorithm in AppendixPBl to check if the system 
is lA feasible. If the network is not feasible, go to Step 6. Otherwise, set d^j = dnj, 
and set V^j*, S^*-, S^* to be matrices aggregated by the basis vectors of V^^, and S^^-, 
SJ^^-, respectively, Vn G {1, G}, j G {1, K}. Exit the algorithm. 

• Step 6 Update stream assignment: Update dg'k' = dgty — 1 and go back to Step 2, 
where {g\ k') is given by (the first line of) (l22l) . 

Remark 4.2 (Subspace Design Criterion in Algorithm^: Similar to the stream assignment 
criteria (|22l) . the core space {V^ } and free space {S* } should be designed to alleviate the 
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lA feasibility constraint as much as possible in order to enhance the network DoF. Hence, 
both (l30l) and (|3TI) are designed to maximize the difference between the number of freedoms 
in intermediate precoder design minus the number of inter-cell lA constraints. ■ 
Remark 4.3 (Relationship between Algorithm^ and\5]): In fact, Algorithm [5] is a back- 
ward compatible extension of Algorithm [2l When the network is fully connected. Step 2~4 
in Algorithm [5] will generate {V^.}, {S^}, and {S^^-} with = S„, Vj G {1,...,^}} 
and rank(S^j) = N^j. However, this particular choice of the core space will not offer any 
additional DoF gain compared to other choices of {V^j} and {S^} satisfying constraint (fT4l) 
in the fully connected case. ■ 

D. Solution of Subproblems 13.21 and 13.31 in Partially Connected Networks 



The solution to Problems 13.21 and 13.31 in the partially connected networks are very similar 
to Algorithm 3 and 4, respectively. Details are omitted to avoid redundance. 

V. Performance Analysis 
A. Symmetric MIMO Cellular Networks with Partial Connectivity 

Definition 5.1 (Symmetric MIMO Cellular Networks with Partial Connectivity): A symmet- 
ric MIMO cellular network has G BSs (each with antennas) serving K MS (each 
with A^^' (< A^*) antennas) per BS. There are at most d^^^"^ = df data streams per MS. 
The partial connectivity is induced by the path loss effects as well as the transmit spatial 
correlation effects due to local scatterings ifTSl . lfT6l . Links from BS-n to MSs of BS-g with 
J < \n — g\ < G — J are assumed to have huge path losses and they are regarded as not 
connected. Furthermore, Ri{< N^) and R2i< N"^) denote the ranks of the intra-cell links 
and inter-cell links. ■ 

As a result, there are three key parameters, i.e. J, Ri and i?2, which characterize the 
connection density, the rank of the intra-cell and inter-cell links of the symmetric MIMO 



cellular network, respectively. In particular, the BS side partial connectivity in Definition 12.11 
of the above symmetric network is given by: 

span ({n(g) : q G Mi(A;)}) , if:g = n 

span ({n(q') : q G 'K2{n — g)}) , if: g ^ n,0 < \n — g\ < J or \n — g\ > G — J (33) 
C^*, otherwise. 



^The transmit spatial correlation is caused by the lack of scattering in the propagation environment around the BSs. 
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where {n(g)}, q G {0, 1, A^* — 1} form a basis for C^', 

Ri{k) = {0, 1, A^* - l}\{A;i?i mod A^*, kRi + 1 mod A^*, {k + l)Ri - 1 mod A^*}, 

M2(m) = {0, 1, A^* - l}\{mi?2 mod A^*, mi?2 + 1 mod N\...,{m + l)i?2 - 1 mod A^*}, 

span({n}) denotes the linear space spanned by the vectors in set {n}. To make sure the 
direct links can have sufficient rank, we also assume: and djK < N*. 

Remark 5.1 (Partial Connectivity in Practice): In practice, singular values of channel ma- 
trices or the path gain of links can hardly be absolutely 0, and hence, the DoF defined by the 
asymptotic slope of the throughput-SNR curve may not correspond to the number of data 
streams transmitted. However, this shall not jeopardize the value of the proposed algorithm, 
i.e. Algorithm [3l HI and [51 This is because in practice, we are interested in the performance 
at finite SNR regime only. As long as the singular values or the path gains are below a 
sufficiently small threshold, we shall quantize the singular values and the path gain to be 
zero and the said channel is partially connected according to Definition 12. 1[ ■ 



B. Analytical Results 

Theorem 5.1 (Achievable DoF of the Proposed Scheme): The total DoF achieved by the 
proposed scheme in a symmetric MIMO cellular network in Definition 15.11 is lower bounded 
by GKd*, where d* is the number of streams assigned to each MS, given by: 

+ ] ] (34) 

mm{G-l,2J)K^ + l'mm{G-l,2J)K + 2j jj ' 

Proof: Please refer to Appendi}<PGl for the proof. ■ 
The following are some interpretations of the results in (|34l) . 

Remark 5.2 (Backward Compatibility with Fully Connected K-pair Interference Channels): 
Consider a special case of fully connected G-pair interference channel with K = 1, J > ^, 



d* = min 


Ri, 


max 1 









G+1 



, which is consistent with 



Ri = R2 = . The achievable DoF in (1341) reduces to 
result in the conventional lA feasibility conditiorO. ■ 
Remark 5.3 (How Partial Connectivity Affects DoF): When the partial connectivity effect 



is strong, i.e. J < i?2 <^ A^*, the network total DoF becomes GK min ( Ri, 



2JX||+1 



Hence, it can be observed that partial connectivity affects the total DoF in three aspects: 



^Using the conventional lA feasibility condition in 1101 for G-pair MIMO interference channels, we have TV* + 
{G+l)d>0^d<i^^. 
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• Gain due to the connection density: As the connection density parameter J limits the 
maximum number of MSs that each BS may interfere, the total DoF of the MIMO 
cellular network is 0{G), which scales with the number of the BS. This behavior 
represents a significant gain compared with the fully connected case in which the total 



• Loss due to weak intra-cell links: Note that the rank of the direct link Ri is one of the 
terms in min function and hence, the partial connectivity may also limit the system DoF 



Remark 5.4 (DoF Scaling Law w.r.t. Number of Antennas): 

• Strong inter-cell link case: When the inter-cell links are strong, i.e. R2 — N^, in (l34l) . the 
second term is the max operation is larger, hence, the total DoF scales on 0{N''' + N*). 

• Weak inter-cell link case: When the inter-cell links are weak, i.e. R2 ^ A^*, in (|34l) . the 
first term is the max operation is larger, hence, the total DoF scales on 0{N'^N^). 

Comparing the two cases, we can see that antennas are more "effective" when the inter-cell 
links are weak. This is because when inter-cell links are weak, the partial connectivity can 
be exploited to eliminate part of the potential interference, thus alleviating the constraints on 
transceiver design. ■ 



In this section, we shall illustrate the performance of the proposed scheme by simulation. 

Definition 6.1 (Randomized Partially Connected MIMO Interference Channels): Consider 
a MEMO cellular network with G BSs and OK MSs. Each BS has A^* antennas and each MS 
has A^** antennas and requests df data streams. The BSs and MSs are distributed uniformly 
in a 30A;m x 30A;m area. All BSs transmit at power P. Denote Dgk^n as the distance between 
the BS n and the A;-th MS of BS-(^. The network is partially connected due to: 

• Path loss effect: If Dgk^n > L, we assume the channel from the BS n to the k-th MS 
of BS-g is not connected, i.e. Hg^ „ = 0. 

• Local scattering effect: If Dgk^n < L, due to local scattering effect, channel fading are 
correlated (only transmit correlation), and hence: 



DoF= C(l) [[TOl. 

• Gain due to weak inter-cell links: When the network is dense, i.e. J ^ 1, 
2LKR2 • Hence, a ^ factor gain can be further observed. 




when the intra-cell links are weak, i.e. small Ri. 



VI. Numerical Results 



gk,n 




gk.n 



{0} 



(35) 
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where Lgfc,„ = span ({eAri(g) : q G {l,...,iV*} and satisfies dSl)}), eAr(cu) = ^[1, e"^'^''^'^) 



1 



< 



sm I 



q_ 

TV* 



< 



A^* - 1 



, G [5gfe,n — -P'a(>S', dgk^n), dgk,n + -P'a(5', (igfc.n)] (36) 



where Fa{S,dgk,n) 



arcsin when: 5" < dg^^r. 



" ' " , S is the local scattering radius 
TV when: S > dgk,n 

as illustrated in Fig. H^, 9gk,n is the angle between the antenna array normal direction 

and the direction from BS n to the k-th MS of BS-g. Please refer to [[T6]| for details. ■ 
The proposed scheme is compared with 5 reference baselines belowOij: 

• Simplified proposed scheme (Baseline (BL) 1): The stream assignment and transceiver 
matrices are designed by Algorithms [21 [3] and HI As we illustrate in Remark 14.31 
Algorithm [2] is a simplified version of Algorithm [51 

• Naive extension of the existing lA algorithm [8] (BL 2): The transceivers are designed 
by naive extension of iterative lA algorithm in [8J as described in Algorithm \T\ 

• Coordinated beamforming |[I7]| (BL 3): The BSs jointly optimize their precoders to 
improve the overall system SINR performance using the algorithm proposed in IfTTll . 

• Round robin scheduling with Intra-cell zero-forcing (BL 4): The BSs are scheduled 
to transmit using round robin. Zero-forcing precoders are adopted. 

• Isotropic transmission (BL 5): The BSs and the MSs apply random precoders and 
decorrelators, respectively. 



A. Fully Connected MIMO Cellular Network 

Fig. [51 illustrates the sum throughput versus SNR (10 log]^Q(P)) for the proposed scheme 
and 5 baselines for an lA feasible MIMO cellular network with G = 3, K = 2, df = 1, 
iV< = 5, A^*" = 2, L, 5" > 30V2km. BL 4 can only achieve 2 DoF as each BS has only | of 
the time to transmit. BL 2 achieves only 3 DoF due to the cross link - direct link overlapping 
issue. The throughput of BL 3 saturates at high SNR since coordinated beamforming [[TTl can 

'The feasible bound for lA algorithms on partially connected MIMO cellular network is still unknown. Therefore, we 
cannot plot the theoretical upper bound as one of the benchmarks in simulation. 

^Note that BL 1 is a simplification of the proposed scheme. It does not address the partial connectivity issue. BL 2 is 
generalized from |8| to the MIMO Cellular Network. Comparison with BLl illustrates the importance of exploiting partial 
coimectivity. On the other hand, comparison with BL2 illustrates the necessity of the decomposition approach proposed in 
this paper. 
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only suppress part of the interference. On the other hand, the proposed algorithm and BL 1 
achieve 6 DoF, which is an achievable upper bound. This result also confirms the comments 
made in Remark 14.31 that Algorithm [5] is a backward compatible extension of Algorithm [2] 
and they have the same performance in fully connected networks. 

B. Partially Connected MIMO Cellular Network 

1) Performance w.r.t. SNR: Fig.[6]illustrates the sum throughput versus SNR (10 log^Q^P)) 
for the proposed scheme and 5 baselines in a MIMO cellular network with G = 12, K = A, 
df = 2, = 8, = A, L = 15km, S = 3km. The throughput of BL 2 also saturates at 
high SNR since the network is not lA feasible. BL 1 achieves 1 1 DoF only as Algorithm [2] 
fails to exploit the benefit of partial connectivity. On the other hand, the proposed algorithm 
achieves 35 DoF, which is significantly higher than all the baselines. Furthermore, the total 
DoF upper bound for the fully connected MIMO network is only 1 1 . This demonstrates that 
partial connectivity can indeed contribute to the significant gain in system throughput. The 
comparison between the proposed scheme and BL 1 illustrates the importance of incorporating 
partial connectivity topology in the lA algorithm. 

2) Performance w.r.t. Partial Connectivity Factors: To better illustrate how different partial 
connectivity factors such as path loss and spatial correlation affect system performance, we 
illustrate the sum throughput versus L (the maximum distance that a BS can interfere a MS) 
and S (the radius of the local scattering) under a fixed SNR (30dB) in Fig. |71 By comparing 
the performance of the proposed scheme with different partial connectivity parameters, we 
have that the performance of the proposed scheme roughly scales O (-^), which illustrates 
a consistent observation as Remark 15.31 that weaker connectivity can indeed contribute to 
higher system performance. Moreover, comparison of the proposed algorithm with BL 1 
further illustrates the importance of adapting the transceiver structures given in Def. 14.11 to 
exploit partial connectivity. By dynamically adapting the transceiver structures, the proposed 
scheme obtains significant performance gain over a wide range of partial connectivity levels. 

Appendices 

A. Proof for Theorem I3.il 

Lemma -A.l (lA Feasibility Conditions of MIMO Cellular Network): With i.i.d. fading. Prob- 
lem 12.11 is equivalent to the following problem almost surely. 
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Problem -A.l (Transformed lA Problem): 



max 



G K 
^ ^ ^nj 



K 



s.t: rank(U;^;t) = dgk, rank([V;i, ...V;;^]) = ^dnj 



(37) 

(38) 
(39) 



e {0, 1, rf--}, Wn,ge{l,...,G},nj^g,k,je{l,..., K}. U 
Proof: We need to show that a) if {dnj,\Jgk,y^nj} is a solution of Problem 12.1 [ there 
must exists {dnj, U^^, V,^j} which is a solution of Problerr FA.il and b) vise versa. 

• Proof of a): Denote the transceivers in the solution of Problem |2. II as {U*^, V*^}. Let 

U'nj = ^nr ^'nj = ^*nr ^^^^ ®' ^ ^^^^ dM]) and ([391). 
. Proof of b) : Denote the solution of Problem PAll as { d;^- , U'^ , V;* }. Note that { U'^ , V;* } 
are functions of the cross link channel states, i.e. {iigk,n '■ 9 7^ n}, which are independent 
of the direct link channel states, i.e. {H„fc „}. Hence we have 



rank 



\ 



L * nl5 •■• * 



/* ] 
riKl 



K 



/ 



Y^dnj, VnG (40) 

i=i 



almost surely. Let U*^ 
from (HOI), Theorem [331 holds, which means {U* , V* } satisfies (E]), ^ and 



U'^, V^* = V;*. and use Algorithm H to construct V;.. Then 



(U^)''H„i,nV„fc = 0, Vn G {1, G}, VjV e {1, K}. 
From ([H]), (I39l) we have: 



(U 



gk) 



H 



,MVn, = 0, Vn ^ G {1, G}, Vj, A; G {1, K}. 



(41) 



(42) 



From ([IB, (ini), (gB, and dH, we have that {U;^-, V;^} satisfy Q^®. ■ 
From Lemmc TATl we need to show A) an optimizing solution {c?*^,, U^*, V^*;^.}, of Probler ri-A.ll 
is a feasible solution of the Subproblems 13. 1[ 13. 2[ 13.31 and B) vise versa. 
First consider the statement A). We first have two lemmas: 

Lemma -A.l (Non-overlapped Subspaces): V is a subspace uniformly distributed in Grass- 
mannian Q{D, N). For any N — D dimensional subspace §, V fl S = {0} almost surely. 

Proof: When < [f J , we have: VnS = {0} ^ v / (§)^, Vv G V ^ ^max(V, (§)^) < 
|, where 6'jnax(''^> denotes the largest principal angle between subspace A and B. Note 
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that both V and (S)"*" are D dimensional subspaces with D < From Theorem 1 in ifTSl . 
we have Pr(e^,.(V, (§)^) < f ) = 1. 

Similarly, when D > [f J, we have: VnS = {0} ^ s / (V)^, Vs G § ^ ^max(§, (V)^) < 
|. Note that both § and (V)"*" are — D dimensional subspaces with N — D < and 
(V)-*- uniformly distributes in Grassmannian Q(N — D, N), we can again apply Theorem 1 
in ifTSl . This completes the proof of the lemma. ■ 

Lemma -A.3 (Uniformly Distributed Precoder Space): Under the i.i.d. fading assumption, 
span(V'*), uniformly distributes in Qi^^^id*^^, N*^, where V* and d*^ are the optimal 
solution of Problem TA-ll 

Proof: In Problen f A.l[ span(V'*) is a function of the channel states {Ylgk^m}, denote 
as span(V';) = F({Hgfc,^}). For any two elements in ^(Ef=i <j> ^n)' V" and V^ denote 
M^. = {{Hgfe,„,} : V^. = F{{Ylgk,m})} = {{Hgfc,„(x)}}, where x G {a, b}. Then there exists 
a unitary matrix T G C^"^^" such that 

yfc _ 'pY'', where is the matrix aggregated by the basis of V^;, x G {a, b}. (43) 

, H„fcm(a.)T, if: m = n, 

Construct a mapping G : W ^ Ir,iigk,rn{b) = \ ■ Substitute 

\Hgk,m{a) , otherwise. 
(1431) into ( |39l ). we have that G is bijective. Denote -D^ and as the probability density 

function of {11^,^ ,„} and V* , respectively. Then from the i.i.d. fading assumption, we have 

that DhiG{{ligk,m{a)})) = Dh{{B.gk,m{a)}). Since G is bijective, we have 

D,{Y')= [ Dh{x)Ax= f Dh{G{x))dx= f Df,{x)dx = D^iY""). (44) 

where x denotes the elements in or H^. With (l44l) . we complete the proof. ■ 
The constraints (l38l) . (I39l) in Problerr TA.ll are the same as that addressed in IfTOl . except 
that the number of transmitter and receiver are different when K > 1. Yet, note that the 
analysis in IfTOll can be easily extended to the case with different number of transmitters 
and receivers, Lemmc TATI enables us to extend the existing lA feasibility conditions to the 
cellular case. Hence from Lemmc TATI and the lA feasibility conditions obtained in [[TOll . we 
have that the feasibility conditions of MIMO cellular network are given by (flOl) . 

Denote {d'^j}, {V*}, {U^\.} as the optimizing variables for Problerr TA.ll Substitute the 
stream assignment policy {d'^j} into Problem |3.U the feasibility condition (flOl) shall be satis- 



fied almost surely. Consider the singular value decomposition of V'^ = A] 



V 

n 





fB 



V\H 
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gk 



are iV* x A^*, Y.^^i <j x E"=i «^nj' ^^fc x ^g'fe. and rf^;-, x unitary matrices, respectively, 



^gk 



A'^ and A'^^ are the first Yl'j=i d-nj and d*f. columns of the corresponding matrices, and 



/U 



[B^k)" 



<.S^^,(BS,)^, where A^, B^, A% B^, 



and S^^, are J2f=i ^nj ^ ZljLi '^nj and d*^ x d*^, diagonal matrices, respectively. Note that 
rank(V*) = J2f=i ^Ij, rank(U*fc) = d*i,, Wn,g, k, we have Sjf and S^;, are full rank. Hence 
we can set: 



A'^ 
A' 



n 

U 
gk 



v:Br(sr) 



(45) 
(46) 



From Lemmc TAJl span (V^*) uniformly distributes in Q{J2j=i ^nj, N^n)- Note that span (V„) 



span (V^*), dim(span (S„)) = iV* — J2j=i ^nj' fro"^ Lemma rA.2[ we have: 

span(V„) n span(S„) = {0}, almost surely. (47) 

From equation (fT4l) . the columns of {V^^} and S„ form a basis for C^". Hence, there 
exist matrices (Ef=i x Ef=i <j) and Q„ ((iV* - Ef=i x Ef=i o^^j) such that: 

Qn 



V = rv^ V 



n2' * niC' 



s„] 



(48) 



From (1471) . R„ is full rank almost surely. Hence we have: 

Set V'f; = Q„R;^^^ where R^^^^ is the matrix aggregated by the (Y.k-i'^njj 



nj 



(49) 



1 to 



(I]i=i ^nj)-th column of R"^ Substitute {V^j}, {11^^} into Problem O and substitute 
{V^^}, {U^;,} into Problem [33l From ©, ([45]), ^ and (|49l), all the constraints are 
satisfied almost surely. This completes the proof for statement A). 

Then to Statement B). From statement A) and [10], under the lA constraint (flOl) . the 
optimal value of Problem 13. 2[ [33] are almost surely. Hence, denote {V*y }, {U*^} 

as the corresponding outputs by solving the three problems sequentially. Substitute these 
outputs to Problem [2J1 then from ([B]), and we have that {d^j}, {V;^}, {U*^} 
must satisfy dH) and ([5]). This completes the proof. 
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B. Low Complexity lA Feasibility Checking 

For notational convenience, denote v^^j, v^^f^ and Cgk^nj as f^, and Cgn, respectively, where 
n = in,j), g = A;), n,ge {1, ...,G}, j,k G 

• Initialize the constraint assignment: Randomly generalize a constraint assignment 
policy, i.e. {c*„g, c^„} such that: c*„g, G NU{0}, 4g + c^„ = Cg^. Calculate {F^, F^}: 

• Update the constraint assignment: As illustrated in Fig. [3l while there exist "over- 
loaded nodes", i.e. P* < or < 0, do the following to update {Cg^, Cgnl- 

- A. Initialization: Select an "overloaded node" with negative pressure. For instance, 
assume < 0, we set P^ to be the root node of the "pressure transfer tree", 
which is a variation of the tree data structure, with its nodes storing the pressures 
at the precoders and decorrelators, its link strengths storing the maximum number 
of constraints that can be reallocated between the parent nodes and the child nodes. 

- B. Add leaf nodes to the pressure transfer tree: 
For every leaf node P^ (x E {t,r}y. 

For every g: If c^g > 0, add P^ as a child node of P^ with link strength c^g, 
where x is the element in {t,r} other than x. 

- C. Transfer pressure from root to leaf nodes: For every leaf node with positive 
pressure, transfer pressure from root to these leafs by updating the constraint assign- 
ment policy {Cg^^, Cg^^}. For instance, as illustrated in Fig. |3^, F*^ "^^S Fg^ ^^"^ 
F*2 is a root-to-leaf branch of the tree (red lines). Update: (c^^gj' = c^^g^ — e, 
(Cgmi)' = Cg,„, + e, (4in2)' = Cg,„2 - (c*„2gj' = ^^gi + e- Hence we have 
{Pk)' = Pk-^ and (F*J' = F*^ + e, where e = min (-F*^ , F*^, ^g^ , c^^„J . 

- D. Remove the "depleted" links and "neutralized" roots: 

* If the strength of a link become after Step C: Separate the subtree rooted from 
the child node of this link from the original pressure transfer tree. 

* If the root of a pressure transfer tree is nonnegative, remove the root and hence the 
subtrees rooted from each child node of the root become new trees. Repeat this 
process until all roots are negative. For each newly generated pressure transfer 
tree, repeat Steps B~D (Please refer to Fig. ISjC! for an example). 

- E. Exit Conditions: Repeat Steps A~D until all trees become empty (hence the 
network is lA feasible) or no new leaf node can be added for any of the non-empty 
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trees in Step B (hence the network is lA infeasible). Exit the algorithm. ■ 

C. Proof for Theorem \3.3\ 

We shall first prove the optimality part via the following two lemmas: 
Lemma -C.l: The updated decorrelators {Ugk} in Step 2 of Algorithm [3] are the optimal 
solution for problems (l24l) . 

Proof: Please refer to ISl for the proof. ■ 
Lemma -C.l: The updated free elements in precoder {V^^} in Step 3 of Algorithm |3] are 
the optimal solution for problems (l23l) . 

Proof: Denote = E^=i,^„ Eti ^".(UffcH...^)^ (Uf,H,fc,„). Note that Q„, is a 
positive semidefinite matrix, and we have in (l23l) : 

G K G K 

trace ((UjlH,fc,„V^^.)''(USH,,,„V^^.)) = J^J^ trace ((V^^O'^Qn.V^.) 

n = l j = l n = l j = l 

#9 ^9 



G 



trace ((q|,V^^. + Qj^.S„V^^.)''(Qiv5 + Q^-S^V 



n = l j=l 
7^9 



= ||QI,VJ,. + Q^^.S„V^,.||| (50) 
By minimizing the Frobenius norm in (l50l) . we have: 

^£ = - ((Qis„)^(Q|^.S„))"' (Qis„)^QtY^^. = -(Sf Q„,S„)-^S^Q„,V^^. (51) 



(pTI) proofs the lemma. 

Now we begin to prove the convergence part. Denote 

G G K K 

^ = E E E E trace ((Uf,H,,,„V^)^(Uf,H,,,„V^)) (52) 

g=l n = l k = l j = l 

#9 

Then / is non-negative. Moreover, from Lemm JCTl I-C.21 / is non-increasing in each 
round of update. Hence, following the analysis in [[8]|, Algorithm [3] is surely to converge. 



D. Proof of Theorem \ 

Lemma -D.l: rank ( [V^t, ^n2^ ^nx] ) = Ef=i dn,. Vn G {1, .., G}. 

Proof: Denote the g-th column in [V^l, V^* , V^}^] and [V^i, V^2, - , ^uk] as v^(g) 
and v,^(g), g G Q = {1, Xl^Li ^nj}, respectively. From the intermediate precoder structure 
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(fT2l) . and the orthogonal constraint (fT4l) (which means that all the v5^(g), g G {1, Z]j=i '^nj} 
and the row vectors in S„ are orthogonal to each other), we have 

vi{q) e span(v^(g)) + span(S„), v^(g) ^ span(S„) (53) 
(span(v^(g)) + span(S,,)) fl {+p^g,eQ (span(v^(p)) + span(Sn))) = span(S„) (54) 

where span(X) denotes the linear space spanned by the column vectors of X. From (|53]) . 
(|54|) . we have y^n{q) ^ span({v/j(p), j> ^ q,E Q}), Vg G Q. This completes the proof. ■ 

Lemma -D.2: Denote the p-th row of U*jH„j „ as w„j(p), p G {1, dnj}- Then we have 
span(w^(p)) follows i.i.d. uniform distribution in Grassmannian ^(1, A'^). 

Proof: Denote h„j_„(g) as the g-th row of iinj,n, Q G {1, A^^^}. From the i.i.d. 
fading assumption, span(h^ „((/)) follows i.i.d. uniform distribution in Q(1,N^). Since in 
(1231) . (1241) . the intra-cell channel states {H„j-„}, n G {l,...,^}, j G {1,...,^} do not 
appear, {H„j_„} are independent of {V,^^} and {U*^}. Hence we have span(w,^(p)) = 
span ( J2q=i '^njip, Q)^nj,nil) ) ^^1^ Uniformly distributed in N^), where Unj{p, q) is the 
element in the p-th row and q-th column of U*y . Note that the rows in U*^ are orthogonal 
to each other, span(w^(p)) are independent p G {1, ...,dnj}- Moreover, since H„j^„, 'Hnj',n 
are independent, span(w;^(j9)) and span(w,^,(j9')) are independent. ■ 

From Lemmc TD^ we can easily deduce the following two results 

K 



rank 



TT_ff TT* TT* "fj^^ TT* 

and if Ylfe=i ^^nj < N^, 

V ^ span ( [ H^i,„U:i, H^s.nU:^, Kk^u^Ik 



^^(inj, almost surely (a.s.) (55) 



Vv G C^"^^ a.s. 



(56) 



If dnj = N^, from Lemmc TDTTI and (1551) . (l25l) is proved. Otherwise, from (l56l) : 



dim ^span 



"trJi TT* "H"^ TT* 

■n-nl,n^nl5 ^nK,n^nK 

K 

Y^dnj + {K-Y,dnj) = Nl a.s. 



(span([V^l,...,V^*^]))- 



dim 

^ AT 



span 



TT-f/ TT* TT^^ TT* 



nspan([V^*i,...,V^*^]) =0, a.s 



H^i,„U:„...,H^.„U;^ ) ^7^([V^^V^*„...,V^*^]) = {O}, a.s. (57) 



where A/^(X) and 7^(X) denotes the null space and the range of matrix X, respectively. From 
(|57]) . we have (1251) . This ends the proof. 
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E. Proof of Theorem 13.51 

We need to show: A) Prove that the output of Algorithm H] {V*^}, n G {1, ...,G}, j E 
{1,...,K}, is indeed a solution for Problem I3.3[ i.e. it satisfies constraints (fT6l) . (flTl) and 
(fTSl) : B) Show that the intra-cell interference power (fTSl) under {V*j} is 0. 

We shall first prove the A) part. In Step 3, from the property of SVD, we have: 

v:f v:^ = (a;^)^a'„^ = I (58) 



a 

'-'rig 



T 



(59) 



where S^^ is a dng x dnq diagonal matrix. 



From Theorem 13.41 in Step 1, L„(q') is full rank almost surely. Note that Rn(q') is a lower 
triangular matrix, in Step 2, we have: 

rank ((U:j^H„,,„V;j = rank(L;(g)) = d„„ G {1, N} a.s. (60) 

{V:X'^np,nKq = 0, ^ 5 G {1, N} (61) 

where L'„(g) is a x matrix aggregated by the elements in the last dng rows and 
columns of L„(g). From (|60l) . we have rank(S^g) = rank(V^g) = dnq almost surely. Hence, 
V*g = V^^^(S^q)~^B„g almost surely. This result leads to the following equations: 

rank ((U:j^H„,,„V:J = rank ((U:,)^H„,,„V:,^(S;j-iB„,) 

= rank ((U:j^H„,,„V;j = Vg G {1, iV} a.s. (62) 

[v:i,v:2,...,v:,,] = [v^t,v^*2,-,v^*;,]-[Qui)(s:,i)-^B„i, 

Q:,(2)(S:,2)-'B„2, QUiV)(S:,^)-^B„^] (63) 

(1581) . (1621) and ([631) proofs the statement A). 

Now we turn to the B) part. Statement B) is an immediate reference of (|6T1) . 



F. Proof for Theorem \4.1\ 

We need to prove statement B) in Appendi^jTAl Following the analysis in ^ Appendix 
B], after introducing the new transceiver structure in Definition 14. 1[ the "no more lA con- 
straints than freedoms" constraint (flOl) is extended to (ITTI) . Moreover, constraints (l28l) . (|29|) 
and the fact that intra-cell channel states {H^^ g} are independent of the inter-cell channel 
states {Hgk.n, n ^ g}, ensures that the statement in Theorem 13.41 still holds under partial 
connectivity. Hence, the solution set of Problem 13.31 is non-empty. Substitute these solutions 
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to Problem [2H then from ([HI), ([B]), and we have that {d*^j}, {Vnj}, {U*^} must 
satisfy (01) and ©. This completes the proof. 

G. Proof for Theorem I5.il 

Due to the symmetry property of the system, in Algorithm [51 the stream assignment policy 
{o?*^}, core space {V^*}, free space {S**^} and linear filter for decorrelator {U^^} hall be 
symmetrical \Jg E {l,...,^}. From dfK < and the first line of (l33l) . in Step 3 of 
Algorithm [5l core spaces assignment is feasible iff. dgk < Ri, j G {1, K}. 

From (1271) . (|3T]) and the second and third lines of (1331) . since the inter-cell partial connec- 
tivity state is the same for all users, in Step 4 and 5 of Algorithm [51 the number of streams 
assigned to each MS dgk, the dimension of the free spaces S"*^,, and the dimension of the 
linear filter Sg^ shall be the same for all MSs. Moreover, as A/^'fc„ = {0}, yn,g,k, we have 
that 5"^;^, shall be the maximum possible value, i.e. — Ri. Denote d = dgk, S = S**^, then 
from the feasibility condition ([TTl) and ([33l) . stream assignment {dgk = d] is feasible if d is 
achievable in the following problem: 

max d (64) 

S.t. min(G - 1, 2J)J^min(d, '^^ ) < ^ + A^'" - d (65) 

d e {0, 1, S e {0, 1, A^* - dK} 
By solving ([641), we get ([341) . 
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Fig. I. Illustration of the cross link - direct link overlapping issue and the partial connectivity in MIMO cellular networks. 
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free space: 




Fig. 2. A simple example of the Core space, Free space and Free elements. In this figure, A'^^ — 2, dgk = 1. 




Fig. 3. Illustrative example of the "pressure transfer tree" and the corresponding operations. A) A tree generated in Step 
A and B; B) Pressure transfer in Step C; C) Removal of depleted links and neutralized roots in Step D. 




Fig. 4. Local scattering effect; The lack of scattering in the propagation environment leads to spatial channel correlation. 
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Throughput Versus SNR under Fully Connected Model (Typology: rf=5,N"=2,G=3,K=2,d =1) 




SNR (dB) 



Fig. 5. Sum throughput versus SNR for the proposed algorithm (and 5 baseUnes) in a fully connected MIMO cellular 
network. 




Fig. 6. Sum throughput versus SNR for the proposed algorithm (and 5 baselines) in a partially connected MIMO cellular 
network. 
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Fig. 7. Sum throughput versus partial connectivity parameters for the proposed algorithm (and 5 baselines). L represents 
the connection density and S is local scattering radius (spatial correlation decreases when S increases). 
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